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Abstract 

Previous studies of connectivity in wireless networks have focused on undirected geometric graphs. More 
^ sophisticated models such as Signal-to-Interference-and-Noise-Ratio (SINR) model, however, usually leads to 

directed graphs. In this paper, we study percolation processes in wireless networks modelled by directed SINR 
graphs. We first investigate interference-free networks, where we define four types of phase transitions and show 
C\| . that they take place at the same time. By coupling the directed SINR graph with two other undirected SINR graphs, 

q ' we further obtain analytical upper and lower bounds on the critical density. Then, we show that with interference, 

percolation in directed SINR graphs depends not only on the density but also on the inverse system processing 
£^ ' gain. We also provide bounds on the critical value of the inverse system processing gain 

I. Introduction 

The study of coverage, connectivity, and capacity in large-scale wireless networks from a percolation 
^ . perspective has attracted much attention recently [l]-[6]. To intuitively understand percolation processes 
in large-scale wireless networks, consider the following example. Suppose a set of nodes are uniformly 
^ . and independently distributed at random over an area. All nodes have the same transmission radius, and 
^ ; two nodes within a transmission radius of each other can communicate directly. At first, the nodes are 
^ , distributed according to a very small density. This results in isolation and no communication among nodes. 
CN ; As the density increases, some clusters in which nodes can communicate with one another directly or 
^ . indirectly (via multi-hop relay) emerge, though the sizes of these clusters are still small compared to the 
^ ; whole network. As the density continues to increase, at some critical point a huge cluster containing a 

•i-H 

^ . large portion of the nodes forms. This phenomenon of a sudden and drastic change in the global structure 
c3 is called a phase transition. The density at which the phase transition takes place is called the critical 
density. A fundamental result of continuum percolation concerns such a phase transition effect whereby 
the macroscopic behavior of the system is very different for densities below and above the critical density 
A c . For A < A c (subcritical), the connected component containing the origin (or any other node) contains 
a finite number of nodes almost surely. For A > A c (supercritical), the connected component containing 
the origin contains an infinite number of nodes with a positive probability [7]— [11]. 

Most previous work employing continuum percolation for the study of large-scale wireless networks 
have focused on geometric models in which a link exists between two nodes when they are within each 
others' transmission radii. This simple model leads to an undirected graph and considerably simplifies 
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the resulting analysis. To reflect realistic conditions in wireless networks, however, more sophisticated 
models for link connectivity can be adopted. For instance, a widely-used model for wireless communication 
channels is the Signal-to-Interference-and-Noise-Ratio (SINR) model [12k [13]. Here, the ability to decode 
the transmitted signal from node i to node j is determined by the SINRL. 

Q.. = W (1 ) 

where P, is the transmission power of node i, is the distance between nodes i and j, and iVo is the 
power of background noise. The parameter 7 is the inverse of system processing gain. It is equal to 1 in a 
narrowband system and smaller than 1 in a broadband (e.g., CDMA) system. The signal attenuation L(dij) 
is a function of distance dij. Under the SINR model, the transmitted signal of node i can be decoded at 
j if and only if /?y > /?, where (3 is some threshold for decoding. In this case, a link is said to exist 
from i to j. Percolation in wireless networks under the SINR model has been studied for the undirected 
case in [4], [5]. Here, it is assumed that the (undirected) link exists if and only if min{/3y, j3ji\ > (3. 
Nevertheless, even if /3y > (3, > (3 may not hold and thus the link (j, i) may not exist. Thus, the 
graph resulting from the SINR model is in general directed. 

Percolation processes in directed lattices have been well studied (see [8] and references therein.) Re- 
cently, based on generating function methods, percolation has been analyzed in directed scale-free random 
graphs [14], and random graphs with given degree distributions [15]. Note, however, that lattices have 
regular geometry, and both scale-free random graphs and random graphs with given degree distributions 
lack the geometric constraints which exist in SINR graphs. Hence the results and analytical methods 
in [8], [14], [15] are not directly applicable for directed SINR graphs. 

To understand how the directional nature of communication links affects the connectivity of wireless 
networks, we first study percolation processes in the SINR model with 7 = (interference-free directed 
graphs), where (directed) link exists if and only if the distance between i and j is less than or 
equal to the transmission radius associated with node i. In such directed graphs, a node has two types of 
links, in-links, which are the links pointing to the node from other nodes, and out-links, which are the 
links pointing out to other nodes. Indeed we can define four types of components with respect to a given 
node u: in-component, out-component, weakly connected component and strongly connected component. 
Corresponding to these four types of components, we can define four types of phase transitions, and further 
define four corresponding critical densities. We will show that all four critical densities are equal to a 
positive and finite value X C (V), which depends on the power distribution at each node. By coupling the 
directed SINR graph with two other types of undirected SINR graphs and using cluster coefficient [16], 
[17] and re-normalization methods [9], we further provide analytical upper and lower bounds on X C {V). 

'Note that the interference term in the conventional definition of SINR is 7 PkL(dkj) rather than 7X^5^ j PkL(dkj) [12], [13]. 

The latter approach assumes that each node knows its own transmitted signal and can subtract it from the received signal. In order to 
keep consistent with previous work [4], [5], we follow this assumption in this paper, though the results of this paper does not rely on this 
assumption. 
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Next, we show that with interference (7 > 0), percolation in directed SINR graphs depends not only on 
the density but also on the inverse system processing gain 7. Indeed there exists a positive and finite critical 
value 7c (A), such that the network is percolated only when A > A C ("P) and 7 > %(A). Furthermore, for 
A sufficiently large, 7 C (A) = O(j). The same results have been obtained in [4], [5] for undirected SINR 
graphs. Our results indicate that the critical inverse system gain has the same asymptotical behavior in 
directed and undirected SINR graphs. 

The remainder of this paper is organized as follows. In Section II, we give the definitions and assump- 
tions for the directed SINR graph model. In Section III, we study directed percolation in wireless networks 
without interference (7 = 0). We formally define four types of phase transitions and their corresponding 
critical densities, and provide analytical lower and upper bounds for these densities. In Section IV, we 
investigate directed percolation in wireless networks with interference (7 > 0). In Section V, we present 
simulation results on directed percolation in SINR graphs, and finally, we conclude in Section VI. 

II. Network Model 

Although some of our results apply to rf-dimensional graphs in general, we will focus on the 2- 
dimensional case. Let || • || be the Euclidean norm, and A = \A\ be the area of A. Assume Xi, X 2 , X„ 
are i.i.d. 2-dimensional random variables with a common uniform density on a 2-dimensional box A = 
[0, ^Jn/X] 2 , where Xj denotes the random location of node i in M 2 . We assume that the transmission 
power Pi are distributed i.i.d. according to a probability distribution fp(p), p G \pmim Pmax], where 
< Pndn < Pmax < °o. This reflects heterogeneity of transmission powers in real wireless networks. We 
further assume 

(i) Pmin > (3N ; and 

(ii) Pr{P = p min } > 0, Pr{P = Pmax } > 0. 

In wireless networks under the SINR model, there is a directed link from node i to node j if > (3, 
and the link is bidirectional if and only if min{/3y, f3ji\ > (3. Denote by G(X n ,V, 7) the ensemble of 
directed graphs induced by the SINR model. In order to show the percolation behavior in G(X n , V, 7), 
we define two other types of undirected SINR graphs: the first is G{X n ,V,^), where there exists an 
undirected link between nodes % and j if and only if min{/3y, PjA- > (3. The second is G'(X n ,V, r )), 
where there exists an undirected link between nodes i and j if and only if m&x{(3ij, f3ji} > (3. The model 
of G{X n ,V,^) was used in [4], [5] as a simplified physical model for wireless communication networks. 

The sum Xlfc^j L(dkj) is a random variable which depends on the locations of all nodes in the network. 
The quantity 

J(x)4 PiL(\\Xi-x\\), xGK 2 . (2) 



SUBMITTED TO IEEE TRANSACTIONS ON INFORMATION THEORY 



4 



is called Poisson shot noise [4], [5], [18]. When is uniformly bounded from below by a nonzero 
constant, the necessary and sufficient condition for J(x) to be finite is 



/"OO 

/ L(x)xdx < oo 

J V 



(3) 



for a sufficiently large y [18]. 

To investigate percolation-based connectivity of directed SINR graphs, we make the following assump- 
tions on the signal attenuation function L(-): 

(i) L(x) < 1 , Vs E (0,oo); 

(ii) L(0) > and 

Prntn 

(iii) L(a;) is continuous and strictly decreasing in x. 

The first assumption reflects the fact that the signal power cannot be amplified by transmitting over 
a wireless channel. With conditions (i)-(iii), ([3]) is guaranteed. Although the last two assumptions are 
introduced for technical convenience, they are also practical in real wireless networks. 

III. Percolation in Wireless Networks without Interference 

For interference-free wireless networks, 7 = 0. To simplify notation, we let G(X n ,V), G(X n ,V) and 
G'(X n ,V) denote G(X n ,V,0), G{X n ,V,0) and G'(X n ,V,0), respectively. 
When 7 = 0, the SINR ® becomes 

Ptj = ^il. (4) 

Thus, in G(X n ,V), there is a directed link from node i to node j if L(dij) > ^0-. Since L(-) is strictly 
decreasing, this condition becomes dy < L^ 1 ^^^j- Let 

R, - L- ^ , (5) 



where r<Ri<r and 

». _ r-l 



iVo/3 



ram 



(6) 



ri'M, (7) 



and 

f^L- 1 ( 

\Pmax . 

Then, G (A^, T 7 ) is the ensemble of graphs where a directed link exists from node i to node j if c?y < 
G(X n ,V) is the ensemble of graphs where an undirected link exists between nodes i and j whenever 
d^ < mm{Ri, Rj}, and G'(X n , V) is the ensemble of graphs where nodes i and j are connected by an 
undirected link whenever dy < max{i?j, Rj}. 

As n and A become large with n/A = X fixed, G(X n ,V), G'(X n ,V) and G(X n ,V) converge in 
distribution to (infinite) graphs G(H\,V), G'(H\,V) and G(H\,V) induced by homogeneous Poisson 
point processes with density A > 0, respectively [10], [11]. According to continuum percolation theory, 
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there exists a positive and finite critical density X C (V) (X' C (V)) for G(H\,V) (G'(H\,V)) such that when 
A > Ac('P) (A > \' C (V)), there is a unique component that contains 9(n) nodes' of G(X n ,V) (G'(X n ,V)) 
a.a.sE This largest component is called the giant component. When A < \ C (P) (A < \' C (V)), there is no 
giant component a.a.s. [9]— [1 1] 

A useful observation is that G G(H\,V) implies G G(H\,V) and G G(H\,V). Also, 
G G(H\,V) implies G G'(H\,V). This relationship will be useful in the following analysis. 

A. Critical Phenomenon 

We now define four types of connected components with respect to a node in G (X n , V). For G (X n , V), 
we use the notation u — > v to mean that there is a directed path from u to v. Similarly, we use the notation 
u <-> v to mean that there is a directed path from u to v and one from v to u as well. 

Definition 1: For a node u G G(Tt\,V), the in-component W in {u) is the set of nodes which can reach 
node u, i.e., 

= {f : f G G(H x ,V),v -> u}. (8) 
The out-component W out (u) is the set of nodes which can be reached from node it, i.e., 

Wout(u) = {v-ve G(H x ,V),u^ v}. (9) 

The weakly connected component W wea k(u) is the set of nodes that are either in the in-component or the 
out-component of node u, i.e., 

W weah {u) 4 W in («) U W out (u). (10) 

The strongly connected component W str ong(u) is the set of nodes that are in both the in-component and 
the out-component of node u, i.e., 

W strong (u) 4 iy m (n) n W out {u). (11) 

Corresponding to these four types of components, there are four types of phase transitions. For instance, 
a directed graph is said to be in the in-component supercritical phase if with probability 1 there exists 
an infinite in-component, and in-component subcritical phase otherwise. In the next subsection we will 
investigate the existence of such a phase transition in directed random geometric graphs. Before that, we 
define the critical densities. Formally, let Ha,o = 7~L\ U {0}, i.e., the union of the origin and the infinite 
homogeneous Poisson point process with density A. Note that in a random geometric graph induced by 
a homogeneous Poisson point process, the choice of the origin can be arbitrary. 

2 We say f(n) — 0(g(n)) if there exists no > and constant Co such that f(n) < cog(n) \/n > no. We say f(n) — Q(g(n)) if 
g(n) = 0(f(n)). Finally, we say /(n) = @(g(n)) if f(n) = 0(g(n)) and /(n) = Q(g(n)). 

3 An event is said to be asymptotic almost sure (abbreviated a.a.s.) if it occurs with a probability converging to 1 as in oo. 
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Definition 2: Let p™(A), p^f (A), p™ eafc (A) and p s £' on9 (\) be the probabilities that the in-component, 
out-component, weakly connected component and strongly connected component containing the origin 
has an infinite number of nodes of the graph G^x^iV), respectively. The critical densities for the 
in-component, out-component, weakly connected component and strongly connected component phase 
transitions are defined respectively as 

A m 4 inf{A:p™(A)>0}, (12) 
\ out 4 inf{A:C(A)>0}, (13) 
\ weak 4 inf{A:p^ fc (A)>0}, (14) 
X strong 4 inf{A:p^(A)>0}. (15) 

By the same argument as in the proof for Theorem 9.19 in [10], it can be shown that when p™(\) > 0, 
G(H\ t o, V) has precisely one infinite in-component with probability 1. It can also be shown that this infinite 
in-component contains a constant fraction of nodes in the network a.a.s. [10]. The same results hold for 
the out-component, weakly connected component and strongly connected component phase transitions. 

Theorem Q] below asserts that all four critical densities of G(H\,V) are actually equal, i.e., all four 
types of phase transitions occur at the same time. The result is not entirely intuitive since it is plausible to 
imagine that a directed random geometric graph experiences three or four steps of phase transitions. That 
is, as the density increases, first an infinite weakly connected component emerges, then an infinite in- 
component and an infinite out-component appear (successively or instantaneously), and finally an infinite 
strongly connected component forms. Nonetheless, we will see that these four types of infinite components 
form at exactly the same time. 

Theorem 1: \ in = \ out = \ wea k — ^strong- 
Proof: Since the choice of the origin in G(H\,V) can be arbitrary, when A > X in , for any node 
u G G(H\,V), Pr{\W in (u)\ = oo} > 0. Because all nodes in G(H\,V) are distributed according to a 
homogeneous Poisson process, for any node u G G(H\,V), the probability Pr{t> G W in (u)} is identical 
for all v G G(H\,V), we have Pr{v G W in {u)} > for all v G G(H X ,V)[_ Equivalent^, for any node 
v G G{H X ,V), Pr{u G W out {v)} > for all u G G{H X ,V), which implies that Pr{| W^O)] = oo} > 0. 
That is A > \ out . Similarly, we can show when A < \ in , then A < X out , and therefore \ in = X out . 

It is obvious that \ wea k < Kn = A out . On the other hand, if A > X wea k, then either p^(A) > or 

> °- Hence we have \ weak = \ in = \ out . 
It is also obvious that X s trong > Kn = ^out- Since the events {|W in (0)| = 00} and {|Wo^(0)| = 00} are 
increasing events @ by the FKG inequality [8]-[10], we have Pr{{|W / i „(0)| = 00} n {1^^(0)1 = 00}} > 

4 Note that, u and v are nodes with fixed labels and random positions. 

5 An event A is called increasing if Ia(G) < Ia(G') whenever graph G is a subgraph of G' , where Ia is the indicator function of A. 
An event A is called decreasing if A c is increasing. For details, please see [8]— [10]. 
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Pr{|W in (0)| = oc}Pr{|W^(0)| = oo}, i.e., p£™»(A) > p£(A)p£*(A). Hence, if A > X in = X out , the 



graph is in the strongly connected component supercritical phase. Thus we have X strong < X in = X out , 



Since all four critical densities are equal, we now define A C (V) as the critical density at which all four 
types of phase transitions take place in G(H\,V), i.e., A C (V) = \ in = X out = X weak = X strong . 

B. Bounds for the Critical Densities 

Instead of directly showing that there exists a positive and finite critical density A C {V) for G (H\, V), 
we provide tight analytical upper and lower bounds on A C ("P). To accomplish this, we couple G (H\, V) 
with G(Hx,V) and G'(H\, V). Due to the relationship between G(TC\, V), G(H X ,V) and G\H\,V), 
it is easy to see that X' C (V) < X C (V) < X C (V). By employing the cluster coefficient method [16], [17] 
and the re-normalization method [9], we proved a lower bound on X' C (V) and an upper bound on X C (V), 
respectively. 

In [16], [17], a new method has been proposed to provide lower bounds on the critical densities for 
d-dimensional random geometric graphs. The methodology is based on the clustering effect in random 
geometric graphs which can be characterized by t-th order cluster coefficients. In the same manner, we 
define the cluster coefficients for G'(H\,V) as follows. 

Definition 3: For any integer t > 3, suppose v 1: . . . ,v t -i £ G'(H\,V) form a single chain, i.e., they 
satisfy the following properties: 

i) For each j = 1, 2, ...,£ — 2, (vj,Vj +1 ) G E. 

ii) For all 1 < j, k < t - 1, (vj, v k ) E for \j - k\ > 1, 

where E denotes the set of links in G'(H\, V). Then the t-th order cluster coefficient C' t is defined to be 
the conditional probability that a node v t is adjacent to at least one of the nodes v 2 , v t -i, given that v t 
is adjacent to v i (averaging over all the possible positions X„ 2 , . . . , X t ^_ 1 in M 2 of the points v 2 , ■ ■ ■ , v t -i 
satisfying conditions (i) and (ii)). 

In general, evaluating C' t , t > 4 for G'(H\,V) is quite difficult. Fortunately, the cluster coefficient C' 3 
can be obtained through a geometrical calculation. The details are given in Appendix A. 

Theorem 2: For 2-dimensional directed SINR graphs, we have for any integer t > 3, 



and therefore A 



strong ^in 



A, 



out- 



□ 



Xc(V)> 



1 



[l-C' t ]vf r g{r)f R {r)dr 



(16) 



where CI is the t-th cluster coefficient for G' 




(17) 
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r'[l - F R (r'))dr', 



(18) 



and 

— > 4 In 2 

V - 6 ^ - 3v ^ ( f- 1} ) r 2 



where <p = sin 1 (~) and r = L 1 ^ 



1 f jM. 

Pm in 



Proof: To show the lower bound on X C (V), we employ the cluster coefficient method [16], [17] to 
obtain lower bounds on X' C (V). Consider a node i E G'(H\, V) at a given location x, with given radius r. 
All the other nodes lying within Afa, r) — the circular region centered at x, with radius r, are adjacent 
to node i. The expected number of nodes in ^4(xj,r) is Xirr 2 . 

A node j within £>(xj,r, f) — the annulus centered at Xj with inner radius of r and outer radius of f, 
is adjacent to % if and only if Rj > \ |xj — Xj| |. The expected number of nodes satisfying this condition is 

/ \2nr' / fn(r")dr"dr' . 

J r J r' 

Thus the mean degree for node i is 

H(r) = Xnr 2 + / X2nr' / fn(r")dr"dr' = \irg(r). 

J r J r' 

The mean degree of G'(H\,V) is 

f/ = Xir / g{r)f R {r)dr. 



Apply the results of Theorem 1 in [16], [17] to G'(H\,V), we have fi' c (V) > jz^- Thus 

A'(^) > w > 

(l-q)7r/;^(r)Mr)rfr' 

which yields the lower bound for A c('P). Because R = L~ x (^§r), we have (fT71) . 

To prove the upper bound, we use a mapping between the continuum percolation model and a discrete 
site percolation model on a triangular lattice to obtain an upper bound on X C (V) and therefore an upper 
bound on A dV). Similar methods have been used in [9]. Let £ r be the triangular lattice with edge length 
r/2. Each site is enclosed by a flower shaped region, which is formed by six arcs of circles. Each of the 
circles has radius of length r/2 and is centered at the midpoint of each edge adjacent to the site. This is 
shown in Figure [TJ 

We say a site is open if there is at least one node of G(TC\,V) in the corresponding flower shaped 
region and closed otherwise. If any two adjacent sites are both open, then the flower shaped regions 
corresponding to these two sites both contain at least one node of G(H\, V). Because R is lower bounded 
by r, the nodes within these flower shaped regions must be directly connected by a link. Consequently, 
if site percolation on the triangular lattice occurs, i.e., there is an infinite cluster of adjacent open sites, 
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then percolation occurs in the continuum model. Since the underlying point process for G(H\,V) is a 
homogeneous Poisson process and flower shaped regions of different sites are disjoint, the probability p Q 
that each site is open is identical and the events are independent of each other. Furthermore, 

Po = l-e- xs f, 

where Sf is the area of the flower shaped region, which can be calculated as 

1 / 3y/3(y/5-l) ^ 2 
Sf = - A [^-^ g Jr, 

where <p = sin -1 (~). 

From the theory of discrete percolation, we know that for site percolation on triangular lattices, the 
critical probability is p c = \ [8]. Thus if 1 — e~ xs f > \, percolation occurs in G(H\,V). Therefore, 

41n2 



K{V) < 



and inequality (1191 ) holds. □ 

A special case of the above model is the directed random geometric graph G(7i\, (a,b)) with binary 
distributed transmission radii, i.e., Pr{_R = a} = p a and Pr{i? = b} = pb, where p a > 0, p& > 0, p a +Pb = 1 
and < a < b. For this model, we calculate the cluster coefficient (detailed analysis is given in Appendix 
B) as 

C = (jp\ + Pi + 3p 2 bPa )C + p b p 2 a ( Ji) " + 1 )(a 2 + b 2 ) + h sin Q x (a + b)} hdh, (20) 

-1 ( h 2 +b 2 -a 2 



where fa = cos 1 ^ fe2 +° 2 fe fc2 j and 9 1 = cos 1 ^ 



2bh 



In light of Theorem [21 we have 



Corollary 3: Let A c (a, b) be the critical density for G(H\, (a, b)), then, 



1 — * 4 In 2 

< A c (a, 6) < t rTiTTT^V- , ( 21 ) 



where <£> = sin -1 (|). 
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IV. Percolation in Wireless Networks with Interference 

Now consider the scenario with 7 > 0. In this case the transmission region of each node is irregular 
instead of circular. Nevertheless, the definitions of the four types of phase transitions are still applicable. 
Percolation in G V, 7) depends not only on A, but also on 7. As before, we can show that the critical 
value of 7 is the same for all four types of phase transitions. That is, in G (H\, V, 7), if 7 is strictly 
less than the critical value, there exists a unique giant strongly (in, out, weakly) connected component. 
Formally, define 

Definition 4: 

7 > C (A) = sup{7 : G(7i\,V^) is percolated, A > A c }, 

7 C (A) = sup{7 : G(Hx,V, / y) is percolated, A > A c }, 

7c(A) = sup{7 : G^HxjV,^) is percolated, A > A' c }, 

where ~\ C {V), \ C {V) and X' C (V) are the critical densities for G(H X ,V), G(H\,V) and G'(H\,V), 
respectively. 

The critical values of 7 depend on the network density A, the distribution of the transmission power 
fp(p), background noise power N , and the threshold (3. To show the existence of "7c(A) and to provide 
lower and upper bounds for it, we employ the same technique as in the previous section. By coupling 
G(H X ,V) with G(H\,V,i) and G'(H\,V,y), it is easy to see that 7c (A) < t" c (A) < j' c (X). Thus, by 
obtaining lower bounds on 7 C (A) and upper bounds on 7c(A), we obtain bounds on 7 C (A). 

Let A c (r) be the critical density for G(H\,r) with r = L~ x ( ). It was shown in [4], [5] that for 

\Pmin J 

any A > A c (r), there exists a constant c\ > such that 7 C (A) > ^. This provides a lower bound for 

tc(A). 

To obtain an upper bound for 7^ (A), we employ a new mapping technique to obtain an upper bound 

on 7^(A). 



Theorem 4: For any A > \' C (V), there exist constants < c 2 < 00 and < c' 2 < A, such that 

c 2 



7 c(A) < 



A 



(22) 



Proof: To prove d22j), we show that 7'(A) < -^h-. Map G'iHx.V^) to a hexagonal lattice C H with 
edge length d > L^ 1 ( ) as shown in Figure [21 Then, L(d) < -2^2.. Therefore, there is no link between 

V Pmax J pmax 

any nodes within two different hexagons that do not share common edges. 

Let the dual triangular lattice be and the hexagon of Ch centered at vertex a of C' T be H a (see 
Figure©. Denote by N(H a ) the number of nodes of G'(TC\,V,j) contained in H a , then N(H a ) has a 

6 The construction of C' T is as follows: let each vertex of C' T be located at the center of a hexagon of Ch- 
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Poisson distribution with mean 



E[N(H a )] = X^d\ 



Now note that if the number of nodes contained in a hexagon is strictly greater than 

Pmax ~ (3N 



N' 



I A 



+ 2, (23) 



y(3p min L(2d) 

then all the nodes in such a hexagon are isolated (i.e, no two nodes share a link). To see this note that 
< L(x) < 1 and L(x) is strictly decreasing in x. If N(H a ) > N', then for any nodes i and j in hexagon 

H a , 

P,L(||X i -X,||) 



Pi. 



^o + 7E*^.,^(l|X*-X 

< 



Pmax 



< 



N + j(N - 2)p min L(2d) 

Pmax 

N + j(N> -2) Pmin L{2d) 
< P, 

and similarly (3ji < (3. This implies that there is no link between nodes % and j. Hence all the nodes in 
the hexagon a are isolated. 

For each vertex a of C' T and hexagon H a of C H , let C a be the event that {N(H a ) < N'}. We call 
hexagon H a and vertex a open if C a occurs, and let 

Po 4 Pr{C a } = Pr{N(H a ) < N'}. (24) 

Let 



d</rf^Km^ 

then 

e 2 E[N(H a )] = 9 2 \^d* - 5A 



(25) 



2 A' C (P)' 
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Fig. 3. An infinite component in G' (7ix,'P,y) implies a path passing through an infinite number of open vertices in £' T (open hexagons 
in C' H ). 



Now choose d sufficiently large so that (1 - 9)E[N(H a )\ - 3 > 0. Let 

Pmax ~ PN 



72 



(3 Pmm L(2d)[(l - 6)E[N(H a )} - 3]" (26) 
Since p max > (3N and (1 - B)E[N{H a )\ - 3 > 0, 72 > 0. Moreover, as (1 - B)E[N(H a )\ > N', by the 

Chebyshev's inequality, when 7 > 72, 

Po = Pr{N(H a ) < N'} 

< Pr{N(H a )<(l-e)E[N(H a )]} 
Var[N(H a )} 



< 



8 2 E[N(H a )} 2 
1 

9 2 E[N{H a )] 

KCP) 

5A 



< - 



where the last inequality is due to A > \' C (V). 

Note that if there is an infinite component in G'(H.x, V, 7), there must exist a path passing through 
an infinite number of open vertices in C T (open hexagons in £#), as illustrated in Figure |3l This is 
because along the infinite component in G'CHxjV,^), each hexagon of Ch contains at least one node of 

Now choose a path in C T starting from the origin having length m. Because the status (i.e., open or 
closed) of each vertex a depends only on the number of nodes contained in the hexagon H a , for different 
vertices a and b, events C a and Cj are independent. Thus 

Pr{30 p (m)} < £(m)p™, 

where O p (m) is a open path in Ct starting from the origin with length m, and £(m) is the number of 
such paths. For a path in C T from the origin, the first edge has six choices for its direction, and all the 
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other edges have five choices for their directions as illustrated in Figure |2] Therefore, we have 

£(m) < 6 • 5 m -\ 

Consequently, 

Pr{3G p (m)} < 6 • 5 m ~V™ = \{$ Po ) m . (27) 

When 7 > 72, p < | and hence |(5p Q ) m converges to as m — > oo. This implies that there is no infinite 
path in C' T a.a.s., and therefore there is no infinite component in G'(J-l\,V,'y) a.a.s. either. By setting 

2V3( Pmax - PN ) 
° 2 9{l-6)f3 Pmin L{2d)<P 1 j 

and 

d ~ ^ (29) 
we complete our proof. □ 

From this theorem, we obtain the following corollary which shows that %(A) and 7 C (A) have the same 
asymptotic behavior with respect to A [4], [5]. 

Corollary 5: For sufficiently large A, 7 > C (A) = {¥). 

V. Simulation Studies 

In this section, we present some simulation results on percolation in directed SINR graphs. Figure |4] and 
Figure \5\ illustrate percolation processes in G(H\, (a,b)) where the transmission powers have a binary 
distribution (see Section III-B). 
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(a) subcritical (b) supercritical 



Fig. 5. Percolation in directed SINR graphs without interference where transmission powers have a binary distribution such that Pr{i?i = 
1} = 0.8 and Pr{i?; = 2} = 0.2: (a) A = 0.75, (b) A = 1. 



Analytical Lower Bound 

4.5 - - Analytical Upper Bound - 

<> Simulation 




1 1.5 2 2.5 3 3.5 4 4.5 5 

b(a=1) 



Fig. 6. Critical density for directed SINR graphs without interference where transmission powers have a binary distribution such that 
Pr{Rt = 1} = 0.5 and Pr{i? 4 = b'} = 0.5 where b' ranges over [1, 5]. 

Figure [6] shows numerical and simulation results for the critical density A c (a, b). Note that the lower 
bound is quite tight. 

In Figure [7J the transmission radius of G(H^\V) has a power law distribution as fn(r) = cr~ a ,r £ 
(1, 2), where a = 3 and c is a normalizing factor. 

In Figure [H simulation results of percolation in directed SINR graphs G(H\,V,'y) with interference 
are shown. In this case, the transmission power at each node has a uniform distribution over [p m in = 
^-iPmax = 2], the background noise power is A^ = 0.1 and the successful decoding threshold is (3 — 0.25. 
The path-loss function is 

L{d - ] = ( d - + ym) 

which satisfies our assumptions on L(-) (i)-(iii). 
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(a) subcritical (b) supercritical 

Fig. 7. Percolation in directed SINR graphs without interference where transmission radii have a power law distribution: fii(r) = cr 
r G (1, 2) with a = 3: (a) A = 0.75, (b) A = 1. 




(a) subcritical (b) supercritical 



Fig. 8. Percolation in directed SINR graphs where the power of each node has a uniform distribution over [p m in = l,Pmax = 2]. 
The background noise power is No = 0.1. The successful decoding threshold is (3 = 0.25. The path-loss function is L(dij) = (dij + 
(2N f3)- 1/:i )- s , and the node density is A = 4: (a) 7 = 0.25, (b) 7 = 0.1. 

In all these figures, a randomly picked source node u is represented by a black solid circle. The nodes 
in Wstrongiu) are represented by red solid circles. The nodes in Wi n (u)\W strong (u) are represented by 
yellow solid circles. The nodes in W out (u)\W str ong(u) axe represented by green solid circles, and those 
nodes not connected to u in anyway are represented by red empty circles. In all of these cases, we see 
that phase transitions with respect to the in-component, out-component, weakly connected component and 
strongly connected component do take place at the same time. 

The values of the critical inverse system gain 7 C (A) depends on the node density and path-loss function. 
Figure [9] shows simulation results for 7*c(A) with node densities ranging from to 5, and with two different 
path-loss functions. In the simulation for G (Ti\, P, 7), P has a uniform distribution over \p m i n , p ma x} with 

Pmin 1 and Pmax 2. 
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0.35 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

Network Density X 



Fig. 9. Critical inverse system gain for percolation in directed SINR graphs as a function of node density, for two path loss functions. 

VI. Conclusions 

In order to understand how the directional nature of wireless communication links affects global 
connectivity, we investigated percolation processes in wireless networks modelled by directed SINR graphs. 
We first studied interference-free networks (7 = 0), where we defined four types of phase transitions 
and showed that they take place at the same time. By coupling the directed SINR graph with two other 
undirected SINR graphs and by using cluster coefficient and re-normalization methods, we further obtained 
analytical upper and lower bounds on the critical density X C (V). Finally, we showed that with interference 
(7 > 0), percolation in directed SINR graphs depends not only on the density but also on the inverse 
system processing gain 7, and there exists a positive and finite critical value 7^ (A), such that the network 
is percolated only if A > X C (V) and 7 > 7c (A). We obtained new upper and lower bounds on 7^ (A). 

Appendix 

A. Calculation of C 3 

Let the transmission radii of nodes i, j and k be r i5 77 and r k respectively, i.e., 

Without loss of generality, we assume that 77 > 77. Denote the coverage area of node i located at Xi 
with radius r« by *A(xj, r^). The cluster coefficient is equal to the conditional probability that nodes i and 
j are adjacent given they are both adjacent to node k. For different ordering on the values of r^, 77 and 
77, the cluster coefficient of G'(7i\,V) is different. In the following, we categorize all the possibilities 
into three scenarios: r k > 77 > r i5 77 > r k > r^, and r, > r, > r k , and calculate the conditional cluster 
coefficient separately. 

Case I: r k > 77 > 
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In this case, to determine the cluster coefficient, assuming both nodes i and j are adjacent to node k, 
the probability that nodes i and j are also adjacent is equal to the probability that two randomly chosen 
points in a circle with radius r k is less than or equal to a distance rj apart. That is, the probability that 
there is a link between nodes i and j is equal to the fraction of A(x.j,Tj) that intersects A(x k ,r k ) over 
A(x k ,r k ). 

This fraction is 

K x ,y) = rin ST 

\A{x k ,r k )\ 

By averaging over all points in A(x k , r k ), we obtain the cluster coefficient for this case: 

1 



C lr k >r 3 >n) = l^^yi / b( Xl y)dxdy 



1 

2 



b(x, y)dxdy. 



Changing to polar coordinates, we obtain 

'{r k >Tj>ri) 



C (r k > rj > ri ) = ^2 I I b{h)hdtpdh. 



k JO JO 

As shown in Figure 10, b(h) can be calculated as 



b[h) = J - — 2 , (30) 



where 

h 2 + r? — ri 



S\ — Zx k XiB = cos 1 



2hrj 



and 

a _ n _i fh 2 + r 2 k -r] 



9i = ZxiX k B = cos 
That is because the area of intersection is given by 



2hr k 



x k BCD XjBAD vx k uxj 
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2 20, 2 20! 1 

= ^. — +^.. — -2.-^81110, 

= 0i^| + Q\t\ — hrj sin 0!. 
Note that b(h) is independent of ip. Hence, 

c {rk> _ r]> _ ri) = -L r r mwdh 

nr k Jo Jo 



r l Jo 



b(h)hdh 
2 f r ' k 

— 4 / (^i r j + ^i r fc ~~ hrj sin <f>i)hdh. 
nr k Jo 



Case II: rj > r k > 

This scenario is similar to the first case, with the roles of r,- and r k exchanged. Thus, we have 

2 f rj 

C (r 3 >r k >n) = — ( d i r l + ^i r i - hr k s m 0±)hdh. 
nr j Jo 

Case III: r 3 - > r,i> r k 

In this case, the probability that there is a link between nodes % and j is equal to the fraction of A(x.j, r,) 
that intersects A(x. k ,ri) over A(x-j,rj). 
This fraction is 

b( x ,y) = — ~rh xi • 

\A*v r i)\ 

By averaging over all points in A{x. j: rj) and changing to polar coordinates, we obtain the cluster 
coefficient for this case as 

c{ r]> _ n> _ rk) = — 2 / mwdh- 

j Jo Jo 

As in Case I, we can calculate 

, , , x 02^ 2 + - hrj sin 6 2 
b (h) = £-2 , 



where 



and 



cf> 2 = Zx jXk B = cos ^ — j 
^ f h 2 + rj - rj 



02 = Zx k XjB = cos 



2/ir,- 



Hence, 



<Wr fc) = i r [ w Kh)hdHh 

2 f/ ' 6(/i)/id/i 



r 2 , 



2 f Tj 

— — I / (02 r j 2 + #2?"f - sm9 2 )hdh. 

*™j JO 
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All of the above results are based on the assumption that rj > r\. Since nodes % and j are equivalent, 
the results also hold when r t > Ty Therefore, the cluster coefficient for G'(H^\ R) can be calculated as 

C {r ] >r k >r l )fR{ r j)fR{ r k)fR{ri)dr j dr k dr i + 
C'ir^n^fR^^fR^dfRir^drjdridrk j . 



r pr rr 



r_ Jn Jr^ 
r pr pr 



r_ Jr k Jri 



B. Calculation of C 

To calculate C, assume both nodes % and j lie within the transmission range of node k. We calculate 
the probability that nodes % and j are also adjacent. Since the roles of nodes % and j are interchangeable, 
we further assume that r,- > 

It is clear that C(r k = a,rj = a, r^ = a) = C(rk = b,r-j — b,ri — b) = C, where C = 1 — is 
the cluster coefficient for random geometric graphs with identical radii [16], [17], [19]. It is also easy to 
check that C(r k = b,rj = b,r { = a) = C(r k = a,rj = 6, = b) = C. Therefore we need to consider 
only two cases: r k = b, rj = a 1 r i = a and r k = a, rj — b,ri — a. 

The first case corresponds to Case I in Appendix A. Thus we have 

C(r k = b,r-j — a,ri — a) — / (0ia 2 + 6\b 2 — hb sin &i)hdh, 

where 



vrfe 2 Jo 

fh 2 + a 2 - b 2 \ 
0i = Zx k XiB = cos -1 ( — J 



and 



,'h 2 + b 2 - n 2 

Ui = /.XiX k B = cos 



2bh 

The latter case corresponds to Case II in Appendix A, thus we have 



2 f b 

C(r k = a,r-j = b,ri = a) = —— / (<fiib 2 + tiia 2 — ha sin tii)hdh. 

™ Jo 

All the above results are based on the assumption that rj > r,. Since nodes i and j are equivalent, the 
results also hold when > rj. Therefore, the cluster coefficient for this model is 

C= (pl+pl + 3p 2 b Pa)C + PbP 2 a (J^j • [(0i + ^i)(a 2 + b 2 ) + hsmti^a + b)] hdh. 
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